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SUMMARY 


Laminated composite material construction is gaining popularity within 
industry as an attractive alternative to metallic designs where high strength at 
reduced weights is of prime consideration. This has necessitated the development 
of an effective analysis capability for the static, dynamic and buckling analyses 
of structural components constructed of layered composites. Theoretical and user 
aspects of layered composite analysis and its incorporation into C.S.A.R. Corpora- 
tion's proprietary version of the NASTRAN* program, CSA/NASTRAN*, are discussed. 
The availability of stress based and strain based failure criteria is described 
which aids the user in reviewing the voluminous output normally produced i ;h 
analyses. Simple strategies to obtain minimum weight designs of composite . uc- 
tures are discussed. Several example problems are presented to demonstrate the 
accuracy and user convenient features of the capability. 


INTRODUCTION 


As structural designers turn more often to high strength light weight 
composite materials to solve critical design problems, the internal loads and 
stress analysis tasks become more complex. The finite element models generated 
to describe composite laminates using existing NASTRAN capabilities either require 
a considerable number of membrane elements ’‘stacked 1 * on top of one another to 
represent the plies or require some form of lumping/ del umping procedure v^en using 
a single element to represent the plies. Both methods of modeling are used to 
overcome the substantial amount of work involved in determining the material 
properties to be referenced by an element. Neither method is without its drawbacks. 
The "stacked" membrane model requires a considerable number of elements and neglects 
bending effects while the single element representation requires pre-processing 
functions to generate an appropriate set of element and material properties and also 
requires post-processing to obtain the individual ply stresses and strains. The 


*CSA/NASTRAN is an advanced proprietary version of the NASTRAN* general 
purpose structural analysis program. NASTRAN is a registered trademark of NASA. 
CSA/NASTRAN is developed, maintained and marketed by the C.S.A.R. Corporation, 
Northridge, California. 
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method of composite material construction analysis discussed In this work and 
Incorporated Into CSA/NASTRAN employs the single element technique. The user 
describes the laminate by means of new PCOMP, PC0MP1, PC0MP2 and MAT8 bulk data 
cards. 

The finite elements used to model the composite structure must be capable of 
representing coupling between the membrant and bending actions to handle the 
general case of unsymmetrlcal ply layups a;id/or element plane offset from the 
plane of the grid points. Since none of the existing NASTRAN combination 
membrane/bending plate elements consider coupling, some work on the finite elements 
available was also required to Implement ths new layered composite anaysls capa- 
bility. After some review of the current e'.ement formulations and Implementations, 
a decision was made to Incorporate new qeneral purpose combined membrane/bendl ng 
plate elements. The two new elements added to the CSA/NASTRAN library are the 
CQUAD4 and CTRIA3. The CQUAD4 Isa 4-node bilinear Isoparametric general quadri- 
lateral element. It Is capable of membrane/bending coupling, variable thickness 
over the element surface and considers the effects of transverse shear flexibility. 
The CTRIA3 Is the 3-node triangular shaped companion to the CQUAD4 element. 

User convenient features to scan the voluminous output normally produced In 
such analyses are provided. The evaluation of a "failure Index" for each element 
based on the commonly used failure theories (maximum stress, maximum strain. Hill, 
Hoffman and Tsal-Wu) allows the user to review at a glance whether any laminate 
is stressed (or strained) beyond allowable limits. Some simple techniques to 
obtain minimum weight designs of layered composite structures are also discussed. 


THEORETICAL DISCUSSION 


Before discussing the theoretical details of the implementation of the layered 
composite analysis capability, the theory of the CQUAD4 and CTRIA3 elements is 
discussed briefly. 


The CQUAD4 E'.ement 

The CQUAD4 element is a four-node bilinear isoparametric element capable of 
representing membrane, bending (with transverse shear effects) and membrane-bending 
coupling behavior. Geometric and kinematic data for the element are shown In 
figure 1. 


Shape of the Element 

Standard Isoparametric theory Is used to represent the shape of the element. 
A iCt of element parametric co-ordinates (;5,"9) have been selected which vary 
linearly between zero and one with the extreme values occurring on the sides of 
the quadrilateral. Lines of constant ^ and lines of constant^ are Indicated on 
figure 1. 
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Element Co-ordinate System 


The x-axIs is along the line connecting the first two grid points; the y-axls 
Is perpendicular to the x-axIs and lies In the "plane" of the element. If all four 
grid points do not lie In a plane, a mean plane Is defined as discussed below. 
Finally, the z-axis Is normal to the plane of the element and forms a right-handed 
coordinate system with the x- and y-axes. 


Mean-Plane 

If the four grid points of the CQUAD4 element do not lie In a plane, a mean- 
plane containing the projections of the four points Is defined such that the four 
points are alternately H units above and H units below the mean-plane as shown In 
figure 2 and reference 1. 


Membrane Behavior 

An enhanced formulation of the four noded Isoparametric quadrilateral membrane 
element available In the MASTRAN^ program, the CQDMEM1 element, is used to represent 
the membrane behavior of the CQUAD4 element. The enhancement consists of using 
reduced Integration to the In-plane shear representation (single point integration, 
at the center of the element, instead of 2x2 integration). All other details of the 
element formulation are the same as that discussed in reference 2. 

If the element Is non-planar, mean-plane transformations that produce only 
forces and not moments at grid points are used to expand the 8x8 stiffness matrix 
to 12x12 to allow for three displacements per grid point. 


Bending Behavior 

A simple, Inexpensive to formulate and accurate bending behavior is a necessary 
prerequisite for plate elements that are to be used in a layered composite analysis. 
The four-node bilinear isoparametric element discussed by Hughes and Tezduyar 
(ref. 3} possesses these qualities and is therefore used to represent the bending 
behavior of the CQUAD4 element. Detailed derivation of the element stiffness, 
load vector calculations and the stress resultants a~e provided in reference 3; 
salient points from the derivations are discussed below. 

The concept is to have the transverse displacement Interpolated via nine-node 
Lagrange shape functions and the »otations via four-node bilinear shape functions. 
The transverse shear strains are calculated In a special way independent of the 
mid side and center node displacement degrees of freedom; hence four-node bilinear 
shape functions may be used for transverse displacement also. In other words, 
the use of special calculations for transverse shear strains allows the use of 
bilinear shape functions for transverse displacement, rotations and transverse 
shear strains together with the benefit to the element of quadratic accuracy with 
respect to Klrchoff modes. 
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The Implementation of the element follows standard isoparametric theory. 
Interested readers may consult references 3 and 4 for additional details. 

If the element is non-planar, mean plane transformations are derived to 
ensure that the element is in equilibrium when the stiffness matrix is transformed 
from grid points on the mean-plane to the actual grid point locations. 


Membrane/Bending Coupling Effects 

The membrane and bending actions are decoupled for plate theory. However, 
practical situations may necessitate use of CQUAD4 elements to model bending of 
plates about an axis offset from the geometrical neutral axis. Membrane-bending 
coupling effects have to be included to analyze such models; these effects also 
occur in unsymmetric laminates. The CQUAD4 element is therefore designed to in- 
clude membrane-bending coupling behavior. 


The CTRIA3 Element 

The CTRIA3 element is a three-node linear element capable of representing 
membrane, bending (with transverse shear effects) and membrane-bending coupling 
behavior. The element and the element co-ordinate system are shown in figure 3. 
The membrane behavior is modeled using the TRMEM element formulation (reference 
2), the bending and membrane-bending coupling behavior are modeled using a pro- 
cedure analogous to that used for the CQUAD4. Being a linear element, the CTRIA3 
element is not as accurate as the bilinear CQUAD4 element. It muse, therefore, 
be mentioned in this context that the CQUAD4 element is to be used for all plate 
modeling requirements; the CTRIA3 element is to be used only where geometrical 
considerations preclude the use of the CQUAD4 element. 


The PSHELL Bulk Data Card 

The properties for both the CQUAD4 and the CTRIA3 elements are specified 
using the new PSHELL element property bulk data card. The PSHELL card data input 
format provide* for specification of element thickness, moment of inertia para- 
meter, transverse shear thickness parameter, stress recovery coefficients and 
material property references. Provisions are made for specifying up to four 
different material property identification numbers to separately represent mem- 
brane, bending, transverse shear and coupled membrane/bending behaviors. 


Force-Displacement Relationship 

The relationship between forces and strains used for the CQUAD4 and CTRIA3 
elements is 
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membrane forces per unit length 


bending moments per unit length 


transverse shear force per unit length 


membrane strains In mean plane 


curvatures 


transverse shear strains 
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and 


G] Is the 3x3 elasticity matrix for membrane action 

G 2 is the 3x3 elasticity matrix for bending action 

G 3 Is the 2x2 elasticity matrix for transverse shear act. 

G 4 Is the 3x3 elasticity matrix fc.' membrane-bendi ng coupling action. 

t Is the element thickness 

I Is the element moment of Inertia 

t s Is the effective thickness for transverse shear 


LAYERED COMPOSITE ANALYSIS 


Composite laminates have a number of laminae stacked at various orientations 
(see figure 4 and reference 5). Ideally, a full 3-dimensional analysis using 
anisotropic material properties Is to be performed. However, lamination theory 
is a good starting point to perform 2 -dlmensional analysis that gives satisfactory 
results at a much reduced cost. 


Assumptions Used in Lamination Theory 
The following assumptions are us r d In lamination t .<;ory: 

(1) Lamina is In a state o r plane stress. 

( 1 i ) Perfect bonding exists between layers so that in slippage of one 
lamina relative to another occurs. 

(Ill) There Is no z-varlation of the transverse displacement; in other words, 
thin plate theory can be used. 

All of the above assumptions are found to be reasonable In practice (espe- 
cially In cases where the thickness of the laminate Is 'mail In comparison with 
the length and width). 

For unidirectional composites, two orthogonal planes of symmetry exist. One 
pla^e Is parallel to the fibers and the other is transverse to the fibers. Know- 
ing the material properties In this system, the 3x3 elasticity matrix [Gr] can be 
evaluated for each lamina. Knowing the orientation of each lamina, the elasticity 
matrix [G^] can be transformed to [Gfr] In a common element system. Knowing the 
elasticity matrix [Gr] for all laminae In the common element coordinate system, 
the membrane, bending membrane-bending and transverse shear elasticity matrices 
for the laminate are calculated as follows: 

[«,] ■ 
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[gJ = Y J 

K] - t* 


The [G3] Matrix is calculated by assuming that the equations of equilibrium 
similar to the simple beam theory can be developed independently for the X- and Y- 
di recti ons. It is to be noted that this is an approximation and that the inter- 
laminar stresses evaluated by lamination theory are only approximate, however, 
it is felt that this approximate analysis is better than an analysis neglecting 
the effects of transverse shear matrix [G3] - 


USER INPUT 


Composite laminate analysis requires input of laminate data and the ortho- 
tropic material property information. This is accomplished by designing the 
following new bulk data cards shown in the Appendix: 

MAT8: two-dimensional orthotropic material data 

PCOMP, PCOMr'l , PC0MP2: Property cards for composites 

User specifies the property identification number of a PCOMP (or PC0MP1 or 
PCuMP2) property card instead of PSHELL property card for the CQUAD4 and CTRIA3 
elements for use as a layered composite element. PCOMP, PC0MP1 and PC0MP2 cards 
refer to the material properties of each lamina using MAT1 , MAT2 01 MAT8 cards. 


OUTPUT FROM LAYERED COMPOSITE AN, .LYSIS 

The stresses and strains output for each lamina for all the elements yield 
voluminous output. Some of these output items are shown In figure 5(a) and 5(b) 
(stresses and strains in each lamina). The concept of failure Index for each 
laminate is Introduced to review the output easily. Five commonly used failure 
theories are provided for this purpose. These are (i) maximum stress theory, 

( 1 i ) maximun strain theory, (iii) Hill's theory, (i«i Hoffman's theory, and 
(v) Tsal-Wu theory. Based on the user specified fa.iure theory, the state of 
stress (or strain) in each lamina is used to evaluate a failure index for the 
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laaina. The hiqhest failure index value aaong all 1 »inae of the laainate and 
the interlaainar stress to allowable bonding stress ratio is defined as the 
failure index of the laainate. By exaaining the failure Index table, the results 
of the analysis can be easily reviewed to see whether any laaina has failed 
according to the spe:1fied criterion. A sample output of the failure Index 
table is shown in figure 5(c). 


NASTRAN MODIFICATIONS 

The design of the layered coaposite analysis capability that has beer in- 
corporated into Ci*/ NASTRAN was driven by two iaportant requireaents. The first 
was that the Cip.bility has aaxiaua versatility and convenience in describing the 
coaposite lauinate with a ainiaua of user action necessary. The second was that 
aodifications to existing NASTRAN functional capabilities be kept to a minimum. 

To achieve the first requireaent, seven new bulk data card types shown in the 
Appendix, and a new functional aodule were designed and added to the program. 

The second requireaent was satisfied without corsiderable effort when the decision 
was aade to iapleaent new finite eleaent technology, rather than u> try and adapt 
the existing memorane/bending plate eleaents to the layered coaposite environment. 
This decision also reduced the risk of inadvertantly disturbing soae aspect of 
the existing finite eleaent iapleaentations. Thus, vhe modificaticis aade to the 
NASTRAN program for tr.e layered composite analysis capability can be divided Into 
two parts: (1) those required for the addition of the new general purpose QUAD4 

and TRIA3 finite eleaents, and (2) those required for the specification and data 
recovery of the coaposite laainate Itself. 

The incorporation of the new 4-node quadrilateral shaped p 1 ate and 3-node 
triangular shaped plate finite eleaents into NASTRAN has been achieved with the 
addition of only three new bulk data cards. These new cards a rt the finite 
eleaent connection cards, CQUAD4 and CTRIA3, and the element property card, PSHELL, 
whi^h is referenced by both elements. This new general purpose capability required 
the following modifications to existing areas of the NASTRAN program for Statics 
and Normal Modes solution sequences: 

° Preface IFSiP and IFXiBD routines to process the CQUAD4, CTRIA3 and 
PSHELL bulk data cards 

° GPTABD block data routine to add the internal descriptions of the new 
finite elements to the NASTRAN element Horary and establish various 
element-dependent pointer data 

° TA1 table processor module routines to process the material property 
orientation ngle definition options available on the C0UAD4 and CTRIA3 
element conu.ction cards 

° EMG module EMGPRO routine to call the new element matrix generation 
subroutines 


74 



° SSG1 module TEMPL routine to call the new element thermal load vector 
generation subroutines 

° SDR2 module routines to call the new element stress data recovery 
subrouti nes 

° OFP module routines to provide additional output file page headings for 
the CQUAD4 and CTRIA3 element forces and stresses 

Modifications to the 0SMG1 module and RANDOM modules are presently underway to 
extend this capability to differential stiffness/buckling and the dynamics 
solution sequences. 

The incorporation of the composite laminate specification capability has 
been achieved with the addition of only four new bulk data cards. These new 
cards are the MAT8 card which describes the material properties of ply layers, 
and the PCOMP, PC0MP1 and PC0W>2 cards which describe the laminate ply layup 
with varying degrees of generality and convenience. One new functional module 
was also required which performs composite element ply stress data recovery 
operations. A preface processor converts the supplied PCWI* and MAT8 data 
into equivalent PSHELL and MAT2 data, updates element connection data 
references to property data and adds the PSHELL data to the element property 
table (EFT) and the MAT2 data to the material property table (MPT) . The preface 
then also sets a DMAP parameter to control data recovery operations based upon 
whether composite element properties are referenced. The layered comoosite 
analysis capability required the following modifications to the NASTRAN code: 

0 Preface IFSiP and IFXiBD routines to process the PCOMP, PC0W>1, PC0MP2 
and KAT8 bulk data cards 

° A new preface processor, IFP6, which generates equivalent PSHELL a:d 
MAT2 bulk data from the PCOMP and MAT8 data supplied; IFP6 also updates 
the GE0M2 data block to reflect the element references to the newly 
generated property data and adds the data to the EPT and MPT 

° XMPLDO routine to establish module properties for the new composite 
e 1 ement stress data recovery processor, SDRCOW 

° FREMAT/MAT material property processing routine to handle the MAT8 data 

3 SDRCOMP, a new functional module added to generate individual ply 
stresses and strains fnr the la/ered composite elements and to generate 
a failure index table of ply data based upon several available failure 
theories 

° OFP module to process the new composite element Failure Index Table, ply 
stress data block and ply strain data block 

The action required on the part of the user to access this new layered composite 
analysis capability is considered to be very minimal. The user must insert the 
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appropriate PCOMP and PAT8 data wnich describe the composite laminate Into the 
bulk data deck. Then If the user wishes to have ply data recovery operations 
performed, a STRAIN= case control request aust be present and the SDRCOHP and 
OFP modules aust be altered Into the rigid format. 

Once the Initial capability verification testing process was satisfactorily 
completed, modifications were made to the NASTRAN rigid format data base to 
provide the necessary CHKPNT/RESTAKT functions associated with the addition of 
the new bulk data cards. Quality assurance testing procedures were then completed 
satisfactorily demonstrating the Incorporation of the layered composite analysis 
capability as a general purpose feature of CSA/NASTRAN . 


CAPABILITY VERIFICATION VESTING 


When any new capability Is added to NASTRAN, a series of tests Is perforaed 
to ensure that the capability has been Incorporated properly. These tests are 
designed to ensure that the new capability performs according to design speci- 
fications end that existing functional ce^blllty has not been adversely 
affected by the new features. The capab" ty verification tests performed to 
ensure satisfactory Implementation of the layered composite analysis feature 
were divided Into two categories. The first series of tests were designed to 
validate t.ie proper Installation of the new CQUAD4 and CTRIA3 finite elements. 

They demonstrated the accuracy of the elements under various geometry, material 
property and loading configurations. Several of the test cases were taken from 
those proposed In reference 6. The remainder of the tests were designed to 
demonstrate Implementation correctness and not necessarily theoretical accuracy. 
The te«t problems considered tc be most Important are now briefly described: 

° Patch test - measures the ability of the element to represent constant 
strain states of deformation. The model geometry and results are shown 
in figure 6. 

° Twisted beam - measures the effect of warping on the element. The model 
geometry and results are shown in figure 1 . 

° Equilibrium tests - each column/row of the elemental stiffness matrix was 
summed about some point (say, node one) of the element to ensure that the 
matrix represented a set of forces that were in equilibrium. This test 
was performed on several configurations which Included non-rectangular 
shapes and element warping. In all instances, the resulting summations 
were computational zeroes. An additional set of tests were also run which 
extracted the free- free mode shapes of the plate using unit masses at each 
degree of freedom. This test resulted in extraction of six distinct rigid 
body modes which further guarantees equilibrium of the stiffness matrix. 

Several other elemental tests were performed to ensure that the elements performed 
according to specifications. These tests are the same as those described in 
reference 6 and the twisted ribbon problem described in reference 3. The results 
from these tests are summarized in table 1. 
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ir. addition to the tests perforaed to validate the addition of new finite 
elements, a series of tests were also perforaed to validate the installation of 
the layered coaposite analysis capability into CSA/NASTRAN. Once again, aany 
sample probleas were run to demonstrate that the capability perforaed up to 
design specifications. These tests included exercising all of the various options 
available on the PCOMP, PC0MP1, and PC0MP2 bulk data cards. Two example probleas 
were run to demonstrate the theoretical correctness of the implementation. These 
problems were taken from reference 7. The geometry for the first problem, a 
static analysis of a simply supported square plate subjected to sinusoidally 
varying pressure load, is presented in figure 8. The geometry for the second 
problem, a modal analysis of a simply supported cylinder, is presented in figure 
9. The modal analysis was accomplished using the cyclic symmetry analysis capa- 
bility within NASTRAN. The results of these example problems are presented in 
table 2. 


EXAMPLE APPLICATION 


To illustrate the application of the layered composite capability in CSA/ 
NASTRAN for design of practical composite structures, the example of an optimum 
design aluminum rectangular plate (ref. 8, 9 and 10) is analyzed using 8-ply 
laminated composite made of graphi te/ epoxy (commercially identified by Type T 300/ 
5208). The geometry, loading and finite element modeling is shown in figure 10(a) 
and 10(b). The material properties and allowable stresses used in the example 
are those specified in reference 5. Starting from the optimally designed aluminum 
plate thickness for the various element groupings, the thickness for the composite 
model is obtained in 2 or 3 iterations using the enhanced fully stressed design 
algorithm implemented in CSA/NASTRAN. The thickness for the layered composite 
plate and the aluminum plate are shown in figure 10(c). The mass for the composite 
model is 0.40275 kg (0.8879 lbs) versus the mass of 1.034 kg (2.28 lbs) reported in 
reference 8 and 1.0796 kg (2.38 lbs) reported in reference 9 and 0.9487 kg (2.0915 
lbs) reported in reference 10. It is to be noted that it is possible to obtain 
other optimum designs for the layered composite model. The design shown in figure 
10(c) is a feasible design (design that satisfies all specified constraints). The 
objective here was to find the weight savings possible between layered composite 
construction and aluminum on a practical example problem subjected to stress and 
displacement constraints. The capabilities provided in CSA/NASTRAN allows the 
engineers the option of reviewing such "optimum" designs tc select the most 
economical designs for their needs. 


WORK IN PROGRESS 


Several additions and extensions to the layered composite analysis capability 
are currently underway, ^or versatility, two additional plate elements are being 
formulated and added to the CSA/NASTRAN finite element library. These elements 
are the 8-node Isoparametric general quadrilateral element CQUA08 and its 6-node 
triangular companion, the CTRIA6. Coth elements will be ab">e to reference com- 
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posite material properties. An option to model the composite elements for membrane 
behavior only (without any reference to ‘.ending, transverse shear or membrane- 
bending coupling) to facilitate modeling the upper and lower covers r .f a wingbox, 
for example, is now available. This option will be extended to include property 
optimization methods to consider the thickness or orientatioi of each individual 
ply as a design variable. The IPTPR1 and IPT2R2 modules are being enhanced to 
incorporate evaluation of design sensitivity coefficients so that structural opti- 
mization using mathematical programming teenniques can be performed in addition 
to fully stress design resizing. 


CONCLUSION 


The theoretical formulations c : layered composite analysis capability for 
addition to the NASTRAN program are presented in this paper. The need to have 
simple, accurate plate elements for layered composite elements is discussed. The 
addition of two such elements, the CQUAD4 and CTRIA3 elements, and the layered 
composite capability into the CSA/NASTRAN program, *n enhanced oroprietary version 
of the April 1984 release of NASTRAN, developed and maintained by C.S.A.R. Corpora- 
tion, is described. Verification problems and example application problems to 
illustrate the useful features available to engineers to analyze and obtain opti- 
mum designs of layered composite structures are also presented. 
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T*b’e 1. Sunary of Test Results for CSA/MASTRAN CQUAD4 Plate Element. 
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♦Regular means that element shape has not been Intentionally distorted 
t Theoretical results are taken from a fine i»esh finite element solution 


Remarks: 1. Letter grades are used to indicate the following error percentages: 


A: Error < 2%; B: 21 < Erro- < 10% C. 1 j% < Error < 20% 

D. 20% < Error £ 50%; F. Error > 50%“ 

2. Letter grades show for problem set numbers 3, 5, 7, 8, 11, 12 and 15 have been 
assigned by averaging the absolute error for each problem in the set. 

3. for details of test problem descriptions 1-15, consult the March 1984 MSC/NASTRAN 
Appllcav.en Mote in the MSC/NASTRAN Application Manual of the MacNeal-Schwendler 
Corporation, los Angeles, CA; for problem 16, refe- to "Finite Elements Based Upon 
Mindlin Plate Theory with particular reference to the Four-node Bilinear 
Isoparametric Element" by T.J.R. Hughes and T.E. Tezduyar, Journal of Applied 
Mechanics, September 1981. 
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TABLE 2. Composite Element Accuracy Test Results 


SQUARE PLATE ~ 3tetic 9 

Center poin t c/m -Elect •on of a. simply supported squsre 
plate aue to 4t sinusoidally Varying pressure load 


Num her of 
Plies 

Type *E 
Laminate 

Center Point 

Oe//ce i ion 

Theory 

NASTRAN 

2 

Un symmetric 

0. 01001*1 

O. 0/0056, 

3 

Symmetric 

0.00035341 

0.00031,245 

4 

Unsymmetnc 

0.0004-21 75 

0. 00042CZ2 


C YL / NOER ~ Cyclic Symmetry Norma! Modes 


Harmon i c 

Longi tuJi.no. ( 

FREQUENCY (rad /sec) 

N umber 

Half oua. i re 

Theory 

NASTRAN 

1 

/ 

IB. Q08 

18. 574 


3 

44.405 

44. 41 4 


5 

Cl. 435 

40.8/0 

2 

/ 

9. 31/ 

<?. 377 


3 

27.2/1 

27. 143 


5 

41.842 

41.535 

3 

/ 

C. 274 

4 .23 0 


3 

/8. 8*5 

/S. 844- 


w 

30.804- 

30.444- 

4 

/ 

4.134 

4.&13 


3 

14.578 

14 . 


5 

24.414 

L_ 

24.404 


81 



























00 

<n 


ACCESS PROGRAM EXAMPLE PROBLEM FEBRUARY 1. 1983 CSA/NASTR*. 10/ 1/84 PAGE IS 

NASA LANGLEY OPTIMIZATION CONFERENCE -LA SCHMIT ET AL 

USE COMPOSITES INSTEAD OF ALUMINUM 


ELEMENT 

ID 

PLY 

ID 

STRESSES IN LA 

STRESSES IN FIBRE AND MATRIX 
NORMAL- 1 NORMAL-2 

Y E R E D 

DIRECTIONS 

SHEAR-12 

COMPOSITE ELEM 
INTER— < 4MINAR STRESSES 
SHEAR- 12 SHEAR -2 Z 

ENTS ICQUAD4) 

PRINCIPAL STRESSES (ZERO SHEAR) 
ANGLE MAJOR MINOR 

MAX 

SHEAR 

101 

1 

-8 

941 10E+02 

-9 

63429E+0? 


60732E+02 

1. 63201E+01 

3 08403E+00 

37 

24 

-7 71930E-02 -1 

10561E >03 

1 66840E+02 

101 

2 

-5 

00003E+03 

-3 

20839E+02 


64373E+02 

2 01973E >01 

2. 49703E+01 

-84 

14 

-4 732228+02 -3. 

O4769E+03 

2 28723E+03 

101 

3 

-4 

71923E+03 

-2 

40933E+02 

2 

7862 4E +02 

2 24036E+01 

3 69024E+01 

86. 

43 

—2. 23 66 IE +02 -4 

736306+03 

2 2 364 IE +03 

101 

4 

-2. 

44336003 

• 4. 

39937E+01 

-2. 

29645001 

2. 23379E+01 

4. 96633E+01 

-89 

43 

-4. 37740E+O1 -2 

44S38E+03 

1 . 20090E+03 

101 

3 

2 

44336E+03 

4 

399376+01 

2 

2964 3E+0, 

2 24036E+01 

3 69024E+0 1 

0 

33 

2 4433BF+03 4. 

37740E+U1 

1 20090E+03 

101 

6 

4 

7 > 923003 

2 

40933E+02 

-2 

7B624E+02 

2. 01973E+01 

2 49703E+01 

-3 

33 

4 7363OE+03 2 

23667E+02 

2 236416+03 

101 

7 

3 

00006E+03 

3 

20P39E+02 

4 

64373E+02 

1 63201E+01 

3. 08402E+00 

3 

86 

3. 04769E+03 4 

732226+02 

2 28723E+03 

101 

8 

8 

941 106+07 

9 

83429! »02 

- 1 . 

60732E+02 

0 0 

-3. 7B227E-06 

-32 

76 

1 103616+03 7 

71930E+02 

1 66B40E+02 

102 

1 

1 

B3268E+03 

-9 

1 3424E+02 

1 

33344E+02 

7 16397E+00 

3. 239036+00 

2 

81 

1 B3933E+03 -9 

22074E+02 

1 38O70E+O3 

102 

2 

-4 

082336+03 

-4 

29399E+02 

-3 

2B776E+02 

8. 76106E+O0 

1. S9068E+01 

-81 

93 

-3. 34397E+02 -4 

1 37336+03 

1 90) ABE +03 

102 

3 

-3 

09688E+O3 

-1 

97373E+02 

3 

17263E+02 

9 71B11E+00 

2 33106E+01 

83 

13 

— I. 70362E+02 -3 

923b9E+03 

1 876/6E+03 

10? 

4 

-2 

37724E+03 

-2 

08939E+01 

-1 

9334BE+01 

9 77636E+00 

3 16419E+01 

-89 

33 

-2. 37371E+01 -2 

37740E+03 

1 17B33E+03 

102 

3 

2 

37724E+03 

2 

08959E+01 

1 

9334B6+01 

9 71B11E+00 

2. 33106E+01 

0 

47 

2 37740E+03 2 

07371 E+01 

1 17833E+03 

102 

6 

3 

89608E+O3 

1 

97373E+02 

-3 

17263E+02 

B 76106E+00 

1 39038E+01 

-4 

87 

3 92389E+03 1 

703626+02 

1 87676E+03 

102 

7 

4 

082356 +C3 

4 

293V9E+02 

3 

287766+02 

7 16397E+00 

3 2390 3E +00 

8 

07 

4 137356+03 3 

343976+02 

1 90148E+03 

102 

0 

-I. 

B3260E+O3 

9 

I3424F+02 

-1 

33344E+02 

0 0 

-2 40970E-06 

-67 

19 

9 22074E+02 -1 

839336+03 

1 3B070E+03 


Figure 5(a): Stress Output for Layered Composite Elements 
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ACCESS PROGRAM EXAMPLE PHOSLEM 

NASA LANOLEY OPTIMIZATION CONFERENCE - L. A SCHMIT ET AL 


FEBRUARY 1, I PBS CBA/NA8TRAN 10/ 1/84 


RASE 


J! 



00 

'J 


USE COMPOSITES INSTEAD OF ALUMINUM 


STRAINS IN LAYERED COMPOSITE ELEMENTS <C0UAD4) 

ELEMENT PLY STRAINS IN FIBRE AND MATRIX DIRECTIONS PRINCIPAL STRAINS (ZERO SHEAR > MAX 


ID 

ID 

NORMAL -1 


NORMAL -2 

SHEAR- 12 

ANGLE 

MAJOR 


MINOR 

SHEAR 

101 

1 

-2 3979P2E-09 

-4. 

4SP33PF-04 

1 . 9442S3E-04 

13 1948 

1 2948P3C-09 

-4. 

890787E-04 

3 4S883SE-04 

101 

2 

-1 B4P977E-04 

-2 

P94044E-04 

-4. 444S43E-04 

-41. 4741 

2 0PP030E-04 

-4, 

P02473E-04 

4 90089 IE -04 

101 

3 

-1 77243PE-04 

-1 

10P744E-04 

2. 4001042 -04 

48 98 3P 

1 89P418E-04 

-4, 

141403E-04 

2 7OO9I0E— 04 

101 

4 

-P 2704B2E-09 

-3 

347PP1E-04 

-8. 808P74E-09 

-74. 8438 

1 7P994IE-04 

-P. 

7 8483 7E -09 

4. P83 1 P'/E -09^fi ;r 

101 

9 

P. 2704S7E-09 

3 

347PS4E-04 

a a08P73E-09 

13. 1948 

P. 784840E-09 

-1 

7P9994E-04 

4 P03IPSE-09g U 

101 

4 

1, 77243PE-04 

1. 

10P744E-04 

-a. 4S0I04E-04 

-41. 4741 

4. 141403E 04 

-1. 

89P4I8E-04 

8. 70091 3E-04O ^ 

r 

101 

7 

1. S4P977E-04 

2 

P94044E-04 

4 444843E-04 

48. 923P 

4. P02471F-0* 

-a, 

0PP030E-04 

4. 90009 IE -04 ^ -r 

101 

S 

2 3979P2E-09 

4. 

4SP340E-04 

-1. 9442S3E-04 

-74 8432 

4. 890780E-04 

-i. 

2948P3E-09 

3 48823PE-04 y. >■ 

10S 

1 

7 P9P403E-09 

-4. 

324P1 IE-04 

1. 30ISS3E-04 

10 0429 

1. 0369I4E-04 

-4 

99944 9E -04 

3 7P048PE-04 "j 

102 

2 

-1. 90P447E-04 

-2. 

43P9S3F-04 

-9. 084337E-04 

-42. 3889 

3 I32P1PE-04 

-7. 

082I4PE-04 

9. I07944E-C4 

102 

3 

-1 443790E-04 

-P. 

09800 IE-09 

3. 09 ISO 2* -04 

47 4119 

1 87P792E-04 

-4. 

24P302E-04 

3. 044927E-04 

102 

4 

-P. 033989C-09 

1 

I370S4E-09 

-1 S9P834E-09 

-7P. P973 

1 . 44444PE-09 

-P. 

344P48E -09 

9 4 1 94PPE-09 

102 

9 

P. 0399P0E-09 

-1. 

1370S7E-09 

1 . B9PS30E-09 

10 0*89 

P 344P91E-09 

-1 

444448E-09 

9 419700E-09 

102 

4 

1 443790E-04 

P 

09S004E-09 

-3 09IS02E-04 

-42 3889 

4 84P308E-04 

-1 

B7P792E-04 

3 044927E-04 

132 

7 

1 90P447E-04 

2 

43P983E-04 

9 084337E-04 

47 4119 

7 088170F-04 

-3 

1 32P 1 PE - 04 

9 I07944E-04 

102 

8 

-7 P9P404E-09 

4. 

324P13C-04 

-1. 30ISS3E-04 

-TP P979 

4 999444E-04 

-1 

02491 4E -04 

3 7P0P70E-04 


Figure 5(b): Strain Output for Layered Composite Elements 
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ACCESS FROOftAM EXAMPLE PROBLEM 

NASA LANOLEV OPTIMIZATION CONFERENCE - L. A. SCHMIT ET AL 
UK COMPOS! Tit INSTEAD OP ALUMINUM 


PCSftUARV !. ITS 5 CSA/NASTSAN 10/ 1/S4 


PAQC 


47 


F A *. L U ft E 
ELEMENT FAILURE PLV 
ID THEORY ID 

101 HILi. I 


INDICES FOR LAYERED COMPOSITE 
PPaFAILURE INDEX FOR PLY FB*FA1LURE INDEX FOR BONDINO 
(DIRECT STRESKS) 1 INTER "LAMINAR STRESSES) 

0 . 0010 


0. 0014 


a 

o 

4 

9 

A 

7 


0. 0039 
0 0010 
0. 0003 
0.0003 
0. 0031 
0. 00*4 


0. 0031 
0 . 0001 
0. 0041 
0. 0031 
0, 0031 
0 0014 


S 0. 0190 


ELEMENTS < C 0 U A D 4 > 
FAILURE INDEX FOR ELEMENT FLAO 

MAX UF FP.FB FOR ALL PLIES 


0 0190 


8 


Figure 5(c): failure Index Table for Layered Composite Elements 
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Membrcne Boundary Conditio** : 

U = 2r - C-kX-*- i>io“ s 

Bending Boundary Condit ions : 

>MO' 3 
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Test Results ; 


Element Output Ci.uantt.ig 

Theory 

MAST RAN 

In- Plane Strains €* , £ y , T 

1.0 x/0"* 

1.0 * 10' 3 

S tresses ^ 

1333. 

1333. 

Shear Stress 

400. 

400. 

Sending Monant pmr ante tengf-h 07 k * My 

1.111*10 ' 7 

l./U'M' 7 

Bendtng Moment par amt length Mag 

0.333* 10' 7 

0.333*10 * 

Surface Stresses <T^ (7^ 

*0.U7 

to. 667 

Surface Shear Stress T»j 

to. ZOO 

*0.200 


Ft jure Co. Element Patch Test Results — CQUAD 4* 
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F.mmd 

End 


Z*\Z QUA 04- Mesh 


Length = IZ. 0 m 

W 'Jih = U *» 

Depth ~ 0.32 w 

Tluist — 10° C root tb tip) 

£1 =* 21 . 0 * 10 * N/m* 21 = 0.22 

A pphed Louis ; Unit -forces mt tip 


Tip Direction Lomd 

Dt+p tenement* Landed Direction j 

Theoretical 

NASTRAN 

Tn- plane Sheer 
Out- erf -pie no Sheer 

0. 005424- 
0. 001754- 

0.005373 

0.001721 


F, 


jure 


7. CQUAD4 Tints teJ Beam Test" Results 
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• 20 QUAD4- Elements 

& ending , Membrane «»<l 
Tronsrrrse Shear Effect* 

• Boundary Coruliiianc 

Edge a.b : 

Ui » 0 
& r •» * 0 

Edge, cd : Simply Supported 

Wr = U« * O 
Qj 3 O 

* Wa.lt thi cknets = 0.5 

• 72 5«y««rtt5 


Figure C f. Simply Supported Composite Cylinder Example 
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6.0 


7 



Composite laminate details: 

8-ply Graphite/Epoxy (T 300/5208) 
(0/+45/ -45/90) Symmetric 


Constraints: 

Element Stress: Hill’s Criteria 

Displacements: _+ .02 at center and corners 
of plate 

(a) Rectangular Plate: Loading and Supports. 
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(Numbers denote element IDs. Numbers In circles 
denote element groups by design variable linking.) 


(b) Finite Element Model. 


Figure 10. Rectangular Plate for Optimum Design 
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APPENDIX 


CSA/NASTRAN Bulk Data Card Descriptions 
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CSA/NASTRAN BULK DATA DECK 


Input Data Card CQUAD4 Quadrilateral Element Connection 


Description: Defines a quadrilateral plate element (QUAD4) of the structural model. This Is an 
Isoparametric membrane-bending element. 


Format and Exaasple : 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CQUAD4 

EID 

PID 

G1 

G2 

G3 

G4 

e 



CQUAD4 

1213 

745 

92 

103 

76 

1001 

32.4 


CQ4 





T1 

T2 

T3 

T4 




+Q4 




.78 



.63 

.60 

.70 





Field 


Contents 


EID Element Identification number (unique Integer > 0) 


PID Identification maaber of a PSHELL or PCOMP property card (Integer > 0 or blank, 

default is EID) 

G1,G2,G3,G4 Grid point identification numbers of connection points (Integers > 0, all unique) 

C Material property orientation specification (Real or blank; or 0 Integer 

^ < 1,000.000). If Real or blank, specifies the material property orientation angle 

In degrees. If Integer, the orientation of the material x-axIs Is along the projec- 
tion onto the plane of the element of the x-axIs of the coordinate system specified 
by the Integer value. The sketch below gives the sign convention for . 


T1,T2,T3,T4 Membrane thickness of element at grid points G1 through G4. (Real >_ 0. or blank, 
not all zero. See Remark 4 for default.) 



x element 
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CSA/NASTRAN BULK DATA DECK 


CQUAD4 (Continued) 


Remarks: 


1. Element Identification nuabers aust be unique with respect to AH other eleaent 
Identification nurters. 

2. Grid points G1 through G4 aust be ordered consecutively around the perlaeter of 
the eleaent. 

3. All the Interior angles aust be less then 180*. 

4. The continuation card Is optional. If It Is not supplied, then T1 through T4 
will be set equal to the value of T on the PSHELL data card. 

5. Stresses are output In the eleaent coordinate systea. 
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CSA/NASTRAN BULK DATA DECK 


Input Data Card CTRIA3 Triangular Element. Connection 


Description: Defines a triangular plate element (TRIA3) of the structural 
model. This Is an Isoparametric membrane* bendl ng element 


Format and Example : 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CTRIA3 

EID 

PID 

Gl 

62 

G3 

e 




CTRIA3 

401 

101 

1001 

1201 

1401 

, 

■a 



CTR3 





T1 

T2 

T3 





♦TR3 



.83 

.81 

.74 






Field Contents 

E1D Element identification number (Unique Integer > 0) 

PIO Identification number of a PSHELL property card (Integer > 0 or blank, default Is 

E1D) 


G1.G2.G3 Grid point Identification numbers of connection point* (Integers > 0, all unique) 

B Material property orientation specification (Real or blank; or 0 < Integer < 

1,000,000). If Real or blank, specifies the material property orientation angle In 
degrees. If Integer, the orientation of the material x-axIs Is along the projection 
on to the plane of the element of the x-axIs of the coordinate system specified by 
the Integer value. The sketch below gives the sign convention for$. 


T1,T2,T3 Membrane thickness of element at grid points Gl, G2, and G3 (Real ^ 0. or blank, not 
all zero. See Remark 2 for default.) 

^element 



38 






















CSA/NASTRAN BULK DATA DECK 


CTRIA3 (Continued) 


Remarks: 1. Element Identification ntanbers oust be unique with respect to all other element 

Identification numbers. 

2. The continuation card Is optional. If It Is not supplied, then T1 through T3 
will be set equal to the value of T on the PSHELL data card. 
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CSA/NASTRAN BULK DATA DECK 


Input Data Card MAT8 Material Property Definition, Fora 8 
Description: Defines the Material property for an orthotropic Material 
Format and Example : 

123456789 10 


MATS 

EID 

El 

E2 

^12 

612 

Gl.Z 

G 2.Z 

RH0 


NAT 8 

101 

28. +6 

.5*6 

0.25 

1.+5 

2. +6 

1.5+6 

0.62 

+A8C 



A1 

1 

A2 1 

TREF 

B 

B 

B 

B 

S 


♦ABC 

24.-6 

1.4-6 ] 

70. 

1.2+4 

1 . 4+4 | 

1 

1 . 3+3 

1.1+3 

1.+3 

+0EF 




f 12 








♦DEF 

1.3-4 










Field Contents 

MID Material ID (1,000,000 > Integer > 01 

El Modulus of elasticity In longitudinal direction (also defined as fibre direction or 

1- direction) (Real i 0.0) 

E2 Modulus of elasticity in transverse direction (also defined as Matrix direction or 

2- di recti on) (Real f 0.0) 

Poisson's ratio {^/^for uniaxial loading in 1 - direction). Note that V ll * € l/€2 

uniaxial loading in 2 - direction is related to 1^2» El. E 2 by the relation 
**12^2 * *^21^1- Real. 

G}; Inplane shear Modulus (Real > 0.0) 

•>1,Z Transverse shear Modulus for shear In 1-Z plane (Real > 0.0 or blank) 

g 2,Z Transverse shear Modulus for shear in 2-Z plane (Real > 0.0 or blank) 

RH0 Mass density (Real) 

A1 Thermal expansion coefficient In 1-dlrection (Real) 

A2 Thermal expansion coefficient In 2-directlon (Real) 

(Continued) 


100 

































csa/nastran bulk data deck 


TREF 

H.*c 

Yt.»c 

s 

P12 

Remarks: 


HAT 8 (Continued) 


Thermal expansion reference temperature (Meal) 

Allowable stresses In tension and conpresslon, respectively. In the longitudinal 
direction. Aequ'red If failure Index Is desired. (Real > 0.0) (Default value for 
X c Is X t > 

Allowable stresses In tension and compression, respectively, In the transverse dlrectlo 
Required If failure Index Is desired. !Real > 0.0) (Default value for v c Is T*) 

Allowable stress for Inplane shear (Real > 0.0) 

Structural damping coefficient (Real) 

Interaction term in the tensor polynomial theory of Tsai-Mu (Real). Required If 
failure Index by Tsai -Mu theory Is desired mid If value of Fjj Is different from 0.0. 


1. If a no &2,z values are not supplied, the In-plane shear modulus will be used. 
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CSA/N^Sl RAN BULK DATA DECK 


Input Data Card PCPP Laye-ed Composite Element Property 

Description : Defines the properties of an n-ply c capos' te material laminate 

Format and Eaagle: 

i 23456789 10 


PCIHP 

PID 

n 

ca 

s 

■g 

TREF 

GE 

OPT 


KPW 




s 

H*FF 




♦ABC 



NIDI 

T1 

6i 

SflTI 

MID2 

T2 

02 

SBUT2 


♦ABC 

10 

0.25 

0. 



YES 





45. 


+0EF 


1 

N1U3 

a 

03 

1 

S*T3 j 

: 

WD4 

T4 j 

04 

SfUT4 


♦DEF 



-45. 




90. 




Field Contents 

PID Property ID (1,000,030 > Integer > 0) 

Z 0 Distance fro* the reference plane to the bottom surface (Real) (Default = 

-1/2 the thickness of the element) 

NSH Non-structural aass per unit area (Real) 

St, Allowable shear stress of the bonding Material (allowable inter- Iwinar shear 

stresr). Required if failure Index is desired. (Real > 0.0) 

F.T, Failure theory to be used to test Aether the element falls or not (BCD). The 

following thejrW »re presently allowed: 

HILL for the Hill theory 
HfFF for the Hoffman theory 
TSAI for the Tsai -Mu theory 
STRESS for maximum stress theory 
STRAIN for maximum strain theory 

TREF Reference temperature (Real) 

GE Damping coefficient (Real) 


(Continued) 
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CSA/NASTRAN BULK DATA DECK 


PCONP (Continued) 


OPT 


HUM 

T1 

•l 


SfcJTi 


Lanl nation generation 09 1 Ion co 4k ( integer _> 0 ) 

Code 

0 all plies are specified 

1 ply layap specified Is syunetrlcal; cnly the plies on one side of the 
ce .rllne are specified 

2 a ‘stack* of «nfcrane-oaly elements will Pe generated for each elenent 
that references this PCONP card 

options nay be contained by adding the specific code values 

Material 10 of the varioes plies. NIDI > 0. N , D2......NIPn > 0. or blank. The 

plies are Identified ty serially mnfeerlng then fro* 1 at the bottoa layer. Tlw 
NIOs aest refer to HAT1 . NAT2, or NATO balk data cards. (Integer > 0) 

Thicknesses of the various piles (Real; TI > 0.0. T2, T3,--> 0.0 o- blank) 

Orientation eagle of the longltadlnal direction of each ply with the naterial 
axis of the elenent. (If the naterial angle on the elenent connection card Is 
G.O. the naterial axis and side 1-2 of the elenent coincide). The plies are to 
be nnOerert serially starting with 1 at the bottoa layer (the bottoa layer Is 
defined as the surface with the largest -Z value in the elenent coordinate 
systen). (Real) 

Stress output required (TES) or not (Ml) for the various ;11es (BCD) (Default * 

NO) 


Renarks: 1. The default under MID2. HID3, .... Is the last defined N*T card, in this case M101; 

so for T2, T3, .... all these thicknesses will be equal to Tl. 

2. At least one of the four values (NlOi. Tl. 0^ . SOUTI) Bust be present for a ply 

to exist. 


3. TREF given on the PCIHP card will be used for all plies cf the elenent; it will 
override values supplied on naterial cards for Individual plies. 

4. IE given on the PCIHP card will be used for the elenent; It will override values 
supplied on naterial cards for Individual plies. 
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CSA/NASTRAN 3ULK DATA DECK 


Input Data Card PCIMP 1 Layered Composite Eleaent Property - Alternate For* 1 

Description: Defines the properties of an n-ply coaposlte aaterlal laminate, all plies being of 
equal thickness. 


Format and Example : 


1 

2 

3 

4 

5 

6 

7 

b 

9 

10 

POMP 

PID 

Zo 

N91 

Sb 

| 

MID 

T 

OPT 


jPCDMP 

891 

m 

18.7 

4. *3 

STRAIN 

10 

0.25 


♦ABC 



®1 



04 

SB 

% 


». 

1 

1 

1 

♦ABC 

f 


♦45. 

-45. 

+90. 

*90. 

-45. 

♦45. 

0.0 

i 


Field Contents 

PID Property ID (1,000,000 > Integer > 0) 

Z 0 Distance fror the reference plane to the bottoa surface (Real) (Default * -1/2 

the thickness of the eleaent) 

NSM Non-structural aass per unit area (Real) 

Sj, Allowable shear stress of the bonding aaterlal (allowable Inter-laainar shear 

stress). Required If failure Index Is desired. (Real > 0.0) 

F.T. Failure theory to be used to test whether the eleaent falls or not (BCD). The 

following theories are presently allowed: 

HILL for the Hill theory 
H#FF for the Hoffman theory 
TSAI for the Tsai -Hu theory 
STRESS for aaxlaua stress theory 
STRAIN for aaxiaua strain theory 

HID Material ID of each ply. The plies are Identified by serially timbering thea 

from I at the bottom layer. The MID must refer to MAT1, MAT?, or MATS bulk 
data sards. (Integer > 0) 

T Thickness of each ply. All plies have the same thickness. (Real > o.O) 


(Continued) 
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CSA/NASTRAN BULK DATA DECK 


PCOMPl (Co rtinued) 


OPT Lawl nation generation option codes (Integer ^0) 

Code 

0 all plies are specified 

1 ply layup specified is syMetrical; jnly the plies on one side of the 
centerline are specified. 

2 a ’stack. " of aewbrane only elewents will be generated for each elewent 
that references this PC IMP card 

Options way be combined by adding the specific code values 

Orientation angle of the longitudinal direction of each ply with the Mterlal 
axis of the elewent. (If the raterlal angle on the elewent connection card is 
0.0, the Mterial axis and side 1-2 of the elewent coincide). The plies are 
to be nuw^ered serially starting with 1 at the bottow layer (the bottow layer 
Is defined as the surface with the largest -2 value In the elewent coordinate 
systew). (Real) 

Rewarts : 1. Stress and strain output will be generated for all plies If a case control 

STRESS request is present. See the PCOHP card description for a wethod of 
reducing the awount of output generated. 
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Input Data Card PC>ff > 2 Layered Composite Element Property - Alternate For* 2 


Description : Defines the properties of an n-ply composite Material laminate 


Format and Example : 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

PCIMP2 

PID 

H> 

NSM 

St, 

D 

Mir 


OPT 


PCINP2 

891 

-1.0 

18.7 

E9 

H0FF 



10 



♦ABC 



T1 

©1 

■a 

*2 

a 

e 3 

T4 

04 


♦ABC 

.25 

g 

.25 

n 


-45. 


♦90. 

♦DEF 



T5 

05 

T6 

*6 



T8 

00 

03 


♦DEF 


♦90. 


-45. 


45. 





Field Contents 

P1D Property ID (1,000,000 > Integer > 0) 

Z 0 Distance from the reference plane to the bottom surface (Real) (Default - 

-1/2 the thickness of tiie eleaent) 

N91 Mon-structural mass per unit area (Real) 

Sb A1 Unable shear stress of the bonding Material (allowable inter-1 aainar shear 

stress). Required if failure Index is desired. (Real > 0.0) 

F.T. Failure theory to be used to test whether the eleaent falls or not (BCD). The 

following theories are presently allowed: 

HILL for the Hill theory 
HfFF for the Hoffaan theory 
TSAI for the Tsai -Mu theory 
STRESS for Maximum stress theory 
STRAIN for aaxIaiN strain theory 

MID Material ID of the various plies. The plies are Identified by serially number- 

ing them from 1 at the bottom layer. The MID must refer to MAT! . MAT2, or 
MAT8 bulk data cards. (Integer > 0) 


(Continued) 
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CSA/NASTRAN BULK DATA DECK 


PC0HP2 (Continued) 


OPT LMination generation option codes (Integer > 0) 

Code 

0 All plies are specified 

1 ply layup specified is syMetrlcal; only the plies on one side 
of the centerline are specified 

2 a “stack* of aeabrane only eleaents will be iterated for each 
element that references this PCOHP card 

Options aay be coablned by adding the specific code values 

Ti Thickness of the various plies (Real; T1 > 0.0, T2, T3,..> .0 or blank) 

Orientation angle of the longitudinal direction of each oly with the Material 
axis of the eleaent. (If the Material angle on the eleaent connection card is 
0.0, the Material axis and side 1-2 of the eleaent coincide). The plies are to 
be neabered serially starting with 1 at the bottoa layer i the bottom layer is 
defined as the surface with the largest -Z value in the eleaent coordinate 
sy steaj. (Real) 


Reaarks: 1. At least one of the two values. f Ti,&j) aust be present for a ply to exist. 

2. Stress and strain outout wi 1 ! be generated for all plies if a cate control 
STRESS reguest is present. See the PCOHP card description for a metnod of 
reducing the Mount of stress output generated. 
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CSA/NASTRAN BULK DATA DECK 



Contents 

Property identification ni*ber (integer > 0) 

Material identification number for the membrane (Integer > 0 or blank) 

Default value for the membrane thickness (Real) 

Material identification nwber for bending (Integer > 0 or blank) 

Bending stiffness parameter (Real or blank, default * 1.0) 

Material Identification nuaber for transverse shear (Integer > 0 or blank), 
must be blank unless MID2 > 0) default * MID2. 

Transverse shear thickness divided by the membrane thickness (Real or blank, 
default * .833333) 

Nonstructural mass per unit area (Reall 

Fiber distances for stress computation. The positive direction is determined b 
the righthand rule and the order in which the grid points are listed on the 
connection card. (Real or blank, see Remark 11 for defaults). 

Material Identification nuaber for membrane-bending coupling (Integer > 0 or 
blank, must be blank unless MIDI > 0 and MID2 > 0, may not equal MIDI or MID2) 

1. All PSHELL property cards must have unique Identification numbers. 

2. The structural mass is computed from the density using the membrane thlckne 
and membrane material properties. 


(Continued) 
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CSA/NASTRAN BULK DATA DtCK 


pshell (Continued) 


3. The results of leaving an HID field blank are: 

MIDI Mo membrane or coupling stiffness 

MID2 Me bending, coupling, or transverse shear stiffness 

MI03 MID2 will be used as default for MID3 

MID4 No bendi ng-membrane coupling 

4. The continuation card Is not required. 

5. The structural damping (for dynamics rigid fonaats) uses the values 
defined for the MIDI Material. 

6. The MID4 field should be left blank If the material properties are synwetrl. 
with respect to the Middle surface cf the shell. If tne element centerline 
Is offset from the plane of the grid points, the MID4 field may be used 

for modeling the offset but It Involves laborious calculations that produce 
physically unrealistic stiffness matrices ("negative terms on factor 
diagonal") If done Incorrectly. 

7. This card Is used in connection with the CTRIA3 and CQUAD4 cards. 

8. For structural problems, PSHELL cards may reference MAT1, MAT2 or HAT8 
material property cards. 

9. If the transverse shear material, HID3, references a MAT 2 data card, the G1 
G12 and G22 values will be used to obtain the symmetric 2*2 [G3] matrix for 
the element. 

10. For heat transfer problems, PSHELL cards may reference MAT4 or MATS materia 
property cards. 

11. The default for 21 Is -T/2, and for Z2 is +T/2, wher - * T is the local plate 
thickness. 
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